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Abstract. 

A formalism is presented to evaluate the Sivers function in constituent quark models. A non- 
relativistic reduction of the scheme is performed and applied to the Isgur-Karl model. The results 
obtained are consistent with a sizable Sivers effect and the signs for the u and d flavor contributions 
turn out to be opposite. The Burkardt Sum Rule is fulfilled to a large extent. After the estimate of 
the QCD evolution of the results from the momentum scale of the model to the experimental one, 
a reasonable agreement with the available data is obtained. A calculation of nuclear effects in the 
pi i' extraction of neutron single spin asymmetries in semi-inclusive deep inelastic scattering off ^He is 

D \ also described. In the kinematics of forth-coming experiments at JLab, it is found that the nuclear 

■S^ ^ • effects arising within an Impulse Approximation approach are under control. 
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°° THE SIVERS FUNCTION IN CONSTITUENT QUARK MODELS 

(N 

C^ ■ The partonic structure of transversely polarized nucleons is one of their less known 

features (for a review, see, e.g., Ref. [1]). The work presented here aims to contribute to 

^ \ the effort of shedding some light on it. 

Semi-inclusive deep inelastic scattering (SIDIS) is one of the proposed processes 
to access the parton distributions (PDs) of transversely polarized hadrons. SIDIS of 
unpolarized electrons off a transversely polarized target shows azimuthal asymmetries, 
c^ I the so called "single spin asymmetries" (SSAs) [2]. The SSAs are due to two physical 

mechanisms, whose contributions can be distinguished O, |4i |5i |6|]. One of them is 
the Collins mechanism, due to parton final state interactions (FSI) in the production 
of a hadron by a transversely polarized quark [2]. The other is the Sivers mechanism 
dVO, producing a term in the SSA which is given by the product of the unpolarized 
fragmentation function with the Sivers PD, describing the number density of unpolarized 
quarks in a transversely polarized target. The Sivers function is a time-reversal odd. 
Transverse Momentum Dependent (TMD) PD; for this reason it was believed to vanish 
due to time reversal invariance. However, this argument was invalidated by a calculation 
in a spectator model [jsl], following the observation of the existence of leading-twist 
Final State Interactions (FSI) [|9D. The current wisdom is that a non-vanishing Sivers 



X 



function is generated by the gauge link in the definition of TMDs flOl, lill, li2L I13|1 . 



whose contribution does not vanish in the light-cone gauge, as happens for the standard 




+ h.c. 



FIGURE 1. The contributions to the Sivers function in the present approach. 



PD functions. Recently, the first data of SIDIS off transversely polarized targets have 
been published, for the proton UM and the deuteron [|l5|l . It has been found that, while 
the Sivers effect is sizable for the proton, it becomes negligible for the deuteron, so 
that apparently the neutron contribution cancels the proton one, showing a strong flavor 
dependence of the mechanism. Different parameterizations of the available SIDIS data 
have been published il6L ll7l. ll8|] . still with large error bars. New data, which will reduce 
the uncertainties on the extracted Sivers function and will help discriminate between 
different theoretical predictions, will be available soon. 

This experimental scenario motivates theoretical estimates of this quantity. Since a 
calculation from first principles in QCD is not yet possible, several model evaluations 
exist, in a quark-diquark model [l8|,[lli[l9|]; in the MIT bag model [201]: in a light-cone 
model f2l']; in a nuclear framework, relevant to proton-proton collisions ['22]. In here, a 
Constituent Quark Model (CQM) calculation will be described [23] . CQM calculations 
of PDs are based on a two steps procedure (see, e.g., Ref. t2M)- First, the matrix element 
of the proper operator is evaluated using the wave functions of the model; then, a low 
momentum scale is ascribed to the model calculation and QCD evolution is used to 
evolve the observable calculated in this low energy scale to the high momentum one, 
where DIS experiments are carried out (see, e.g., Ref t2M)- Such procedure has proven 
successful in describing the gross features of PDs (see, e.g., [|25i|26]) and GPDs (see, 
C-g- [1271 12 SIP, by using different CQMs. Similar expectations motivate the present study 
of the Sivers function. 



The formalism of Ref. [[23|1 to calculate the valence quark contribution to the Sivers 
function from any CQM is summarized here. A difference in the calculation of TMDs, 
with respect to calculations of PDs and GPDs, is that the leading twist contribution to 
the FSI has to be evaluated. All the technical steps of our procedure can be found in 
Ref. 12311 . where a workable formula is derived. The main approximations have been: 
i) only the valence quark sector is investigated; ii) the leading twist FSI are taken into 
account at leading order, i.e. only one-gluon-exchange (OGE) FSI has been evaluated 
in the treatment of the gauge link (see Fig. 1); iii) the resulting interaction has been 
obtained through a non-relativistic (NR) reduction of the relevant operator, according to 
the philosophy of constituent quark models [29], leading to a potential V^r. In terms of 
the latter, the Sivers function for the quark of flavor ^, f^^x.kj), can be written (cf. 



Fig. 1 for the labels of the momenta and helicities): 

f2 n ^2; 



ftf{xM) = Zl-ig^^ j dhdh^^5{k+ -xP+)5ChT+qT-kT)J^A\) 



where g is the strong coupling constant, M the proton mass, and 



j^u{d) 



Y^ ^IfS-^l [k3,m3;ki,mi; P -k3 -kunin 



1±T3(3) 



VNR{ki,k3,q) 



X ^sf,s,=- 1 (4 + q, m'j, ■,ki-q,m[;P-k3-ki, m,A . (2) 

Eq. (dl), with ^"^^' given by Eq. (O, provides us with a suitable formula to evaluate 
the Sivers function, once the spin-flavor wave function of the proton in momentum 
space, i.e. the quantity ^sf, is available in a given constituent quark model. A few 
remarks are in order. First of all, the helicity conserving part of the global interaction 
does not contribute to the Sivers function. Besides, one should notice that, in an extreme 
NR limit, the Sivers function would turn out to be identically zero. In our approach, 
it is precisely the interference of the small and large components in the four-spinors 
of the free quark states which leads to a non-vanishing Sivers function, even from the 
component with Z = of the target wave function. Effectively, these interference terms 
in the interaction are the ones that, in other approaches, arise due to the wave function 
(see, e.g., the MIT bag model calculation lllOn '). 

There are many good reasons to use the Isgur-Karl model pO] to test the performance 
of the approach. First of all, the IK is the typical CQM, successful in reproducing the 
low-energy properties of the nucleon. In particular, in the IK model, one expects small 
corrections from terms O [k^ /nP') . Besides, one of the features of the IK model is that 
the OGE mechanism [12911 . which reduces the degeneracy of the spectrum, is taken into 
account. It is therefore natural to study our formalism, based on OGE FSI, within the 
IK framework. Concerning PDs, it has been shown that the IK model can describe their 
gross features, once QCD evolution of the proper matrix elements of the corresponding 
twist-2 operators is performed from the scale of the model to the experimental one 



[|24l 1251 . |26D. Reasonable predictions of GPDs have also been obtained [27], and this 
makes particularly interesting the evaluation of the Sivers function in the IK model. 
The final expressions of the Sivers function in the IK model are rather involved and not 
presented here, since they can be found in Ref. [23]. 

To evaluate numerically Eq. ([I]), the strong coupling constant g, and therefore (J.s{Q^), 
has to be fixed. Here, the prescription introduced in the past for calculations of PDFs in 
quark models (see, i.e., Ref. Il24ll ) will be used. It consists in fixing the momentum scale 
of the model, the so-called hadronic scale /Iq, according to the amount of momentum 
carried by the valence quarks in the model. In the approach under scrutiny, only valence 
quarks contribute. Assuming that all the gluons and sea pairs in the proton are produced 
perturbatively according to NLO evolution equations, in order to have ~ 55% of the 
momentum carried by the valence quarks at a scale of 0.34 GeV^, as in typical low- 
energy parameterizations, one finds, that jUq ^ 0. 1 GeV^ if A^^^ — 0.24 GeV. This 
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FIGURE 2. Left (right): the quantity /{j^ '''{x), Eq. ©, for the u (d) flavor The dashed curve is the 
result of the present approach at the hadronic scale j^q. The full curve represents the evolved distribution 
after standard NLO evolution (see text). The patterned area represents the 1 — (7 range of the best fit of the 
HERMES data proposed in Ref. ilT 



yields as{lJ.Q) / (A-n) ^ 0.13 [24]. For an easy presentation, the quantity which is usually 
shown for the results of calculations or for data of the Sivers function is its first moment, 
defined as follows : 
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d'kT^ftf{xM 



(3) 



The results of the present approach for the moments Eq. ([3]) are given by the dashed 
curves in Fig. 2 for the flavors u and d. They are compared with a parameterization of 
the HERMES data, corresponding to an experimental scale of 2^ = 2.5 GeV^ [|17l1 . The 
patterned area represents the 1 — C7 range of the best fit proposed in Ref. [170. Clearly, a 
different sign for the u and d flavor is found. 

Let us see now how the results of the calculation compare with the Burkardt sum rule 
pil) . which follows from general principles and must be satisfied at any scale. If the 



proton is polarized in the positive y direction, in our case, where only valence quarks are 
present, the Burkardt sum rule reads: 






S 



)=0 



(4) 



where 



S\ 



{kf) 



dx 



dkrpt/ixM) 



(5) 



Within our scheme, at the scale of the model, it is found {k") = 10.85 MeV, (k^) = 
— 11.25 MeV and, in order to have an estimate of the quality of the agreement of our 



results with the sum rule, we define the ratio 

{14) - (ki) ' '■"' 

obtaining r ~ 0.02, so that we can say that our calculation fulfills the Burkardt sum rule 
to a precision of a few percent. 

The magnitude of the results is close to that of the data, although they have a different 
shape: the maximum (minimum) is predicted at larger values of x. One should anyway 
realize that one step of the analysis is still missing: the scale of the model, /ip , is much 
lower than the one of the data, which is Q^ = 2.5 GeV^. For a proper comparison, 
the QCD evolution from the model scale to the experimental one would be necessary. 
Unfortunately, the Sivers function is a TMD PDs and the evolution of this class of 
functions is, to a large extent, unknown. In order to have an indication of the effect of the 
evolution, we perform a NLO evolution of the model results assuming, for the moments 
of the Sivers function, the ones defined in Eq. Q, the same anomalous dimensions 
of the unpolarized PDFs. As described in the previous section, the parameters of the 
evolution have been fixed in order to have a fraction ~ 0.55 of the momentum carried 
by the valence quarks at 0.34 GeV^, as in typical parameterizations of PDFs, starting 
from a scale of /Iq ^ 0.1 GeV^ with only valence quarks. The final result is given by 
the full curve in Fig. 2 for the u and d flavor. As it is clearly seen, the agreement with 
data improves dramatically and their trend is reasonably reproduced at least for x > 0.2. 
Of course a word of caution is in order: the performed evolution is not really correct. 
In any case, an indication of two very important things is obtained: i) The evolution 
of the model result is necessary to estimate the quantities at the momentum scale of 
experiments, as it happens for standard PDs 1241. |25L 1261]: ii) after evolution, the present 



calculation could be consistent with data, at least with the present ones, still affected 
by large statistical and systematic errors. One should notice that the agreement which is 
found is better than that found in other model calculations il9Ll20|l . especially for what 
concerns the fulfillment of the Burkardt Sum Rule. 



THE SIVERS FUNCTION FROM NEUTRON (^HE) TARGETS 

As we have discussed in the previous section, the experimental scenario which arises 



from the analysis of SIDIS off transversely polarized proton and deuteron targets 1141 11511 
is rather puzzling. The data show a strong, unexpected flavor dependence in the az- 
imuthal distribution of the produced pions. With the aim at extracting the neutron in- 
formation to shed some light on the problem, a measurement of SIDIS off transversely 
polarized -^He has been addressed [32], and two experiments, planned to measure az- 
imuthal asymmetries in the production of leading n^ from transversely polarized ^He, 
are forth-coming at JLab [33]. Here, a recent, realistic analysis of SIDIS off transversely 



polarized ^He [13411 is described. The formal expressions of the Collins and Sivers con- 
tributions to the azimuthal Single Spin Asymmetry (SSA) for the production of leading 
pions have been derived, in impulse approximation (lA), including also the initial trans- 
verse momentum of the struck quark. The final equations are rather involved and they 



are not reported here. They can be found in [34]. The same quantities have been then 
evaluated in the kinematics of the planned JLab experiments. Wave functions [35J ob- 
tained within the AVI 8 interaction [36] have been used for a realistic description of 



the nuclear dynamics, using overlap integrals evaluated in Ref. [(370, and the nucleon 
structure has been described by proper parameterizations of data or suitable model cal- 
culations [16, 38]. The crucial issue of extracting the neutron information from ^He data 
will be now thoroughly discussed. As a matter of facts, a model independent procedure, 
based on the realistic evaluation of the proton and neutron polarizations in ^He 139(1 . 
called respectively pp and p„ in the following, is widely used in inclusive DIS to take 
into account effectively the momentum and energy distributions of the bound nucleons 
in ^He. It is found that the same extraction technique can be applied also in the kine- 
matics of the proposed experiments, although fragmentation functions, not only parton 
distributions, are involved, as it can be seen in Figs. 1 and 2. In these figures, the free 
neutron asymmetry used as a model in the calculation, given by a full line, is compared 
with two other quantities. One is: 

where / stands for "Collins" or "Sivers", Aj^^'' is the result of the full calculation, 
simulating data, and J„ is the neutron dilution factor. The latter quantity is defined as 
follows, for a neutron n (proton p) in -^He: 

and, depending on the standard parton distributions, /^'^(x), and fragmentation func- 
tions, Z)^''' (z), is experimentally known (see [34] for details). A^ is given by the dotted 
curve in the figures. The third curve, the dashed one, is given by 

a;; ^ -\ (a^^"' - ippdpA^^pA , (9) 

i.e. ^He is treated as a nucleus where the effects of its complicated spin structure, leading 
to a depolarization of the bound neutron, together with the ones of Fermi motion and 
binding, can be taken care of by parameterizing the nucleon effective polarizations, Pp 
and pn- One should realize that Eq. ((Tj) is the relation which should hold between the 
^He and the neutron SSAs if there were no nuclear effects, i.e. the ^He nucleus were a 
system of free nucleons in a pure S wave. In fact, Eq. ((Tj) can be obtained from Eq. ((9() by 
imposing p„ = 1 and pp = 0. It is clear from the figures that the difference between the 
full and dotted curves, showing the amount of nuclear effects, is sizable, being around 
10 - 15 % for any experimentally relevant x and z, while the difference between the 
dashed and full curves reduces drastically to a few percent, showing that the extraction 
scheme Eq. Q takes safely into account the spin structure of ^He, together with Fermi 
motion and binding effects. This important result is due to the peculiar kinematics of 
the JLab experiments, which helps in two ways. First of all, to favor pions from current 
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FIGURE 3. Left (right) The model neutron Collins (Sivers) asymmetry for n^ production (full) in 
JLab kinematics, and the one extracted from the full calculation taking into account the proton effective 
polarization (dashed), or neglecting it (dotted). The results are shown for z=0.45 and Q^ = 2.2 GeV^, 
typical values in the kinematics of the JLab experiments. 



fragmentation, z has been chosen in the range 0.45 < z < 0.6, which means that only 
high-energy pions are observed. Secondly, the pions are detected in a narrow cone 
around the direction of the momentum transfer. As it is explained in |34|l . this makes 
nuclear effects in the fragmentation functions rather small. The leading nuclear effects 
are then the ones affecting the parton distributions, already found in inclusive DIS, 
and can be taken into account in the usual way, i.e., using Eq. Q for the extraction 
of the neutron information. In the figures, one should not take the shape and size of 
the asymmetries too seriously, being the obtained quantities strongly dependent on the 
models chosen for the unknown distributions [38]. One should instead consider the 
difference between the curves, a model independent feature which is the most relevant 
outcome of the present investigation. The main conclusion is that Eq. Q will be a 
valuable tool for the data analysis of the experiments JBBl] . 

The evaluation of possible effects beyond lA, such as final state interactions, and the 
inclusion in the scheme of more realistic models of the nucleon structure, able to predict 
reasonable figures for the experiments, are in progress. 
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